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Abstract 

We consider multi-hop multiple access (MAC) and broadcast channels (BC) where commu- 
nication takes place with the assistance of relays that amplify and forward (AF) their received 
signals. For a two hop parallel AF relay MAC, assuming a sum power constraint across all 
relays we characterize optimal relay amplification factors and the resulting capacity regions. We 
find the sum rate capacity and the maximum rate for each user in closed form and express the 
optimal rate pair (i?i, R2) that maximizes HiRi -l-/i2-R2 as the solution of a pair of simultaneous 
equations. We find that the parallel AF relay MAC with total transmit power of the two users 
Pi+ P2 ^ P and total relay power Pr is the dual of the parallel AF relay BC where the MAC 
source nodes become the BC destination nodes, the MAC destination node becomes the BC 
source node, the dual BC source transmit power is Pr and the total transmit power of the AF 
relays is P. The duality means that the capacity region of the AF relay MAC with a sum power 
constraint P on the transmitters is the same as that of the dual BC. The duality relationship 
is found to be useful in characterizing the capacity region of the AF relay BC as the union of 
MAC capacity regions. The duality extends to distributed relays with multiple antennas and 
more than 2 hops as well. For the M hop AF relay MAC where the to*'* hop consists of 
relays (possibly with multiple antennas) with total transmit power P„i and source power is Pqi 
the dual BC is obtained by reversing the direction of communication so that the transmit power 
of the BC source is Pm and the relays participating in the i*'' hop have transmit power Pm-i- 
Thus, even though the transmitters and receivers are switched, each link in the original MAC 
is associated with the same transmit power as in the dual BC. 



1 Introduction 



Duality relationships arc an important part of network information theory for the insights that 
they provide and the novel techniques that they inspire. In particular, wireless networks have 
the potential for many elegant duality relationships due to the inherent topological similarities 
between various cooperative communication scenarios. A case in point is the duality between 
the multiple input multiple output (MIMO) MAC and the MIMO BC [1,2]. It is shown that 
a MIMO MAC with a sum power constraint P across all transmitters has the same capacity 
region as its reciprocal network, i.e. the network obtained by switching the transmitters and 
receivers. The reciprocal network for the MIMO MAC is the MIMO BC with the same channels 
and transmit power P available to the BC transmitter. Thus the MAC-BC duality relationship 
unifies the capacity results of the multiple access and broadcast channels. The MAC-BC duality is 
also instrumental in proving the optimality of dirty paper coding on the MIMO BC with Gaussian 
codebooks [2,3], and has played a key role in one of the rare characterizations of the capacity 
region of a non-degraded broadcast channel [4,5]. Moreover, from a computational perspective this 
duality relationship enables the transformation of non-convex rate optimization problems on the 
broadcast channel into convex rate optimization problems on the dual multiple access channel [6] . 
These simplifications lead to efficient resource allocation algorithms on the broadcast channel. 

The MAC-BC duality is especially interesting as it hints at a fundamental reciprocity property 
in a multiuser framework. The reciprocity of a wireless channel is a rudimentary property for point 
to point (single user) communications: i.e., the channel capacity is unchanged if the roles of the 
transmitter and receiver are switched. With multiple nodes the implications of switching the roles 
of transmitters and receivers are not easily characterized. The MAC-BC duality shows that for the 
single hop Gaussian MIMO MAC and BC such a reciprocity holds. It also confirms the intuition 
that the topological isomorphism between reciprocal wireless channels should translate into simple 
relationships between their capacity regions. 

The MIMO MAC and BC model the uplink and the downlink for one-hop communication over 
cellular networks or a similar Wi-Fi network topology. While one-hop communication continues to 
be an active area of research, there is also increasing interest in communication over multiple hops, 
i.e. relay networks, spurred by the apparent remarkable potential of cooperative communications. 
For example, using fixed relays in cellular systems can significantly improve the communication 
capabilities of the network [7]. Various relay strategies have been studied in literature, primarily 
for point to point communications. These strategies include amplify-and- forward [8,9], where 
the relay sends a scaled version of its received signal to the destination; demodulate-and-forward 
(DemF) [9] , where the relay demodulates individual symbols and retransmits those symbols without 
regard to an outer code; decode-and-forward (DecF) [10] where the relay decodes the entire message. 
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re-encodes it and sends it to the destination; compress-and-forward (CF) [11,12] where the relay 
sends a quantized version of its received signal; and estimate and forward (EF) where the relay 
sends a soft estimate of its received symbol to the destination. Of these, amplify and forward (AF) 
is perhaps the most interesting from a practical standpoint as it requires the relays to only scale 
their received symbols. Thus the complexity and cost of relaying, always an issue in designing 
cooperative networks, is minimal for AF relay networks. In addition to its simplicity AF is known 
to be the optimal relay strategy in many interesting cases [13-17]. It is also shown that the power 
efficiency of AF relaying increases with the number of relays in a random network [18]. 

AF relay optimization for dual hop communications has been the focus of much recent research. 
[19-21] consider the case of orthogonal relay transmissions. While orthogonal relay schemes are 
attractive for wideband communications, Marie and Yates [22] have shown that for amplify and 
forward relays, shared bandwidth transmission schemes can provide higher capacity. Marie and 
Yates also find closed form solutions for the relay amplification factor and the point to point 
AF relay channel capacity with shared band transmission. Optimum and near optimum power 
allocation schemes for single branch multi-hop relay networks have been considered in [23]. 

While much of the work on relay networks has focused on point to point communications, 
multiuser relay networks are increasingly gaining attention as seminal work in this area [7,12,24-26] 
has shown the remarkable advantages of multiuser relaying. In [27] gains for AF relays in a multiuser 
parallel network are determined to achieve a joint minimization of the MMSE of all the source signals 
at the destination. Tang et. al. [7] consider a MIMO relay broadcast channel, where a multiple 
antenna transmitter sends data to multiple users via a relay with multiple antennas over two hops. 
They find different algorithms for computing the transmit prccoder, relay linear processing matrix 
and the sum rate under the assumption of zero-forcing dirty paper coding and Gaussian signals. 
Capacity bounds are used to establish that the performance loss is not significant. Capacity with 
cooperative relays has been explored for the multicast problem by Marie and Yates [28,29], for the 
broadcast problem by Liang and Veeravalli [26], and for the mixed multiple access and broadcast 
problem by Host-Madscn [30]. Marie and Yates explore an accumulative multicast strategy where 
nodes collect energy from previoTis transmissions, while Liang and Veeravalli [26] and Host-Madscn 
[30] address the general question of optimal relay functionality which may not be an amplify and 
forward scheme. 

In this work, we pursue two related objectives. The first is to investigate the capacity optimal 
relay amplification factors for two hop multiple access and broadcast channels where communication 
takes place via parallel AF relay links and no direct link exists between source (s) and destination(s) . 
With a sum power constraint on all the relays, we characterize the optimal relay amplification factor 
for the parallel AF relay MAC for all points on the boundary of the MAC capacity region. We 
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obtain the sum capacity and the individual user capacities in closed form and present a couple 
of simultaneous equations whose numerical solution yields the optimal rate pair (i?i,i?2) that 
maximizes fiiRi + /U2i?2 for any /xi,/i2 > 0. The capacity region of the parallel AF relay BC 
optimized over all relay amplification factors is found to be non-convex and is harder to obtain 
directly. However, we identify an interesting duality relationship that allows us to compute the 
capacity region of the parallel AF relay BC as the union of the relay MAC capacity regions. 

The second objective of this work is to identify duality relationships in AF relay networks. We 
obtain a general duality result for multi-hop multiple access and broadcast channels where each 
hop may consist of parallel AF relays and the relays may be equipped with multiple antennas. For 
the two hop ease, wc show that the multiple access channel with total transmit power of all users 
equal to P and total relay transmit power is the dual of the BC obtained when the destination 
becomes the transmitter and the transmitters become the receivers, and the powers are switched 
as well, i.e. in the dual BC, the transmit power is P^ and the total relay power is P. In general, 
for multi hop AF relay networks with parallel AF relays at each hop, and with possibly multiple 
antennas at the relays, the duality result holds when each link is associated with the same transmit 
power in the original MAC and the dual BC channels. 

2 Single User AF Parallel Relay Channel 

We start with the simple illustrative example of a single user (point to point) parallel AF relay 
channel. 

2.1 Point to Point Channel Model: PTP(F, P, G, {D}, Pr) 
The input-output equations for this case are as follows: 

R = Fx + Nr (1) 
y = TV[GD (Fx + Nr)] + n (2) 

where R, F, G, D and are R x R diagonal matrices with the z*'* principal diagonal terms 
fij fit 9i-i di, n-i^R respectively representing the received signal at the z*'* relay node, the channel 
coefficient from the source to the i*'* relay node, the channel coefficient from the z*'* relay node 
to the destination, the relay amplification factor for the z*'* relay node, and the additive white 
Gaussian noise (AWGN) component at the z*'* relay receiver, modeled as an i.i.d. zero mean unit 
variance Gaussian random variable, x is the source transmitted symbol and n is the zero mean 
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unit variance AWGN at the destination node. The power constraints are specified as follows: 



The point to point channel under these assumptions is denoted as PTP(F, P, G, {D}, Pr). We 
use the notation {D} to indicate all feasible relay amplification matrices while we use D to indicate 
a specific choice of the relay amplification matrix. For example, PTP(F, P, G, D, Pr) refers to the 
point to point channel with a given D matrix. 

We assume fixed channels so that all the channel knowledge is available to all nodes. Note that 
while in general individual power constraints at relay nodes may be a more reasonable assumption 
for distributed nodes, we assume a sum power constraint across all relay nodes instead. This 
assumption may be more relevant for fixed wireless networks with slowly varying channels where 
the channel information overhead is small. The sum power constraint is interesting for two reaons. 
First, it provides useful insights into relay optmization. Optimal resource allocation, e.g., optimizing 
relay amplification factors is a challenging problem especially when the relays serve multiple users 
simultaneously, as in multiple access or broadcast channels. The second reason is that the sum 
power constraint enables powerful duality relationships that are also of fundamental interest. 

Finally, throughout this paper we consider real signals, noise, channel and relay amplification 
factors. However, most of the results apply to the complex case as well. For example, for the point 
to point channel, it is easily seen that even with complex quantities, the relay must choose the 
phase of its amplification factor di to cancel the total phase of fi and Qi. Therefore, there is no loss 
of generality in the assumption that fi,gi,di are real. The exceptions will be pointed out in the 
conclusions section. 

2.2 Capacity and Optimal Relay Amplification 

The capacity of the point to point parallel AF relay channel and the optimal relay amplification 
factors have been obtained previously by Marie and Yates [22]. We re- visit the result in order 
to introduce our notation and to illustrate the key ideas that are later applied to multiple access 
channels. Given any choice of relay amplification vector D that satisfies the relay transmit power 
constraint, the resulting channel is an AWGN channel whose capacity is simply 



Therefore, the capacity optimization problem for the point to point channel may be represented as: 



Source Power Constraint: 



Relay Power Constraint: 



E[x'] = P, 

E[||D (Fa; + Nr) = TV(D2(I + PF^)) = Pr. 



(3) 
(4) 





C^^^(F,P,G,{D},Pfl) 



maxC^^^(F,P,G,D,P,j) 



(6) 
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where P is the set of feasible relay amplification matrices: 

V = {D : Tr(D2(/ + PF2)) = P/j} 
The following theorem presents the closed form capacity expression. 
Theorem 1 [Marie, Yates 2004] The capacity of the point to point channel described above is 



log 



R f2 2 

Jk9k 



The optimum relay amplification vector for the point to point relay channel is given by 

D = 7FG (l + PF[i]2 ^ p^G^\ (7) 



where ^ is a constant necessary to satisfy the relay transmit power constraint, and may be expressed 
explicitly as: 



7 



Pr 



\\ Tr (f2G2 (/ + PF2 + PkG2)-2 ^ pp2) j 

The proof of this result is provided by Marie and Yates in [22]. 



(8) 



2.3 Reciprocity of the Point to Point AF Relay Channel 

From the result of Theorem ^ notice that 

C^^^(F,P,G,{D},Pk) = C^^^(G,Pk,F,{D},P) (9) 

The capacity is unchanged if we switch variables as follows: 

P ^ Pr (10) 

Pr ^ P (11) 

F ^ G (12) 

G ^ F (13) 

In other words, the capacity is the same if we switch the transmitter and the receiver as long as 
the powers of the transmitter and the relay are also switched. Fig. ^ shows the dual channels that 
have the same capacity. 

Further, note that the optimal relay amplification vector was found as 

D = 7FG (I + PF^ + PrG^) (14) 
6 




Note that except for the constant 7 this matrix is also unchanged when apphed to either the original 
channel or the dual channel. However, 7 must be different on the two channels because the power 
constraints are different. For the original channel 7 is given by Theorem^ whereas on the dual 
channel we will have 



Tdual 



^ Tr (f^G^ (I + PF2 + PrG2)-2 (/ + PG2) 



P 

(15) 



The following theorem shows that this duality property is actually much stronger as it holds not 
only for the optimal relay amplification matrix D but rather for every feasible D. 

Theorem 2 

C^^^(F,P,G,D,P^) = C^^^(G,P«,F,kD,P) (16) 

Given any relay amplification matrix D for a point to point parallel AF relay channel PTP{F, P, G, D, Pr) 
there exists a dual point to point parallel AF relay channel PTP(G, Pr, F, kD, P) that has the same 
capacity as PTP{F, P, G, D, Pr) and where k is chosen to satisfy the relay power constraint of the 
dual channel. 

Proof: We start with the capacity and relay power constraint expressions of the original and 
the dual point to point channels. 

C^^^iF, P, G, D, Pn) = log (1 + fj^gg^X;' ) , TriB\l + PF^)) = Pr (17) 

C^^^iG, Pn, F, .D, P) = log (1 + ^^:^gggCr; ) , Tr(.2D2(I + PrG')) = P (18) 
Substituting from the power constraint into the capacity expression we have: 

\ Ti(D^(/+PF^)) + Tr(D2G2) J 

( [TV(GDF)]^PPh ^ 

= Tr(D2(I + PF2 + P^G2) 
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/ [Tr(GDF)]^ PP^ ^ 

= |^' + Tr(D2(I + PF2 + P^G2)j (22) 

Thus, C^^^(F,P,G,D,P^) = C^^^(G,P^j,F,AcD,P). ■ 
The proof of Theorem |2l ihustrates an important point. By normahzing the noise power to 
unity and also incorporating the relay power constraint into the normalized channel definition 
PTP(F, P, G, D, Pr) may be represented as: 



Tr(GDF)^ 

y = — X + n. (23) 

^Tr(D2(I + PF2 + PKG2)) 

Comparing this to the normalized representation of the dual PTP(G, Pr, F, kD, P), 

Tr (GDF) JP]^ , , 

y = ^ ' =x + n. (24) 

v/Tr(D2(I + PF2 + PRG2)) ^ ' 

notice that the normalized dual channel is identical to the original channel. 

The appealing feature of this normalized form is that it is unaffected by any scaling of the relay 
amplification matrix D by any constant k. While there is only a unique scaling factor for any 
diagonal matrix such that it will satisfy the relay transmit power constraint, that scaling factor 
does not affect the normalized channel as the power constraint is already accommodated into this 
form. Further, the reciprocity of the point to point channel is also evident from the normalized 
form. 



3 Parallel AF Relay MAC and BC 




Figure 2: Dual Parallel AF Relay MAC and BC channels 

We consider two hop multiple access and broadcast channels where communication takes place 
with the assistance of multiple parallel AF relays. As in the point to point case, all channels are 
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known and fixed, there is one sum power constraint across all relays, and all channel, input, output, 
noise variables and the relay scaling factors are real. For simplicity we focus on the two user case. 
The AF relay MAC and its dual BC models are depicted in Fig. [5J 



3.1 Symbol Definitions and Notation 



Variable 


AF MAC 


Dual AF BC 


R 


AT 1 J? n 1 A T — " 1 

IN umber of parallel At relays 


JN umber of paraflef Al^ refays 


Pi 


Transmit power for source 1 




P2 


iransmit power for source z 




F 


Pi + P2 


iotal transmit power of H relays 


Pr 


iotai transmit power of H relays 


J. "J. 

Source transmit power 


A3] 

J i 


Channef from j source to i relay 


Channel from « relay to j destmation. 


9i 


Channef from t relay to destmation 


Channel from source to i relay. 


di 


Amplification factor at the i^^ relay 


Amplification factor at the i^^ relay 


X 




Transmitted symbol for common source 


Xi 


Transmitted symbol for source i 




y 


Received symbol at common destination 




Vi 




Received symbol at i^^ destination 




Received symbol at the i^^ relay 


Received symbol at the i*^ relay 


diri 


Symbol transmitted by the i*^ relay 


Symbol transmitted by the i*'' relay 




Unit power AWGN at the i^^ relay 


Unit power AWGN at the i*^ relay 


m 




Unit power AWGN at i*'^ destination 


n 


Unit power AWGN at common destination 





3.2 Two user Parallel AF Relay MAC and BC Models 

We use the notation MAC(fW, Pi, pl^l , P2, G, {D}, Pr) to denote the two user multiple access 
channel described above, i.e. with the channels between the transmitters and relays f[^],F[^], the 
corresponding source transmit powers Pi, P2 (respectively), the channel between the relays and the 
destination G, the relay amplification vector D, and total transmit power at all relays Pr. To 
distinguish the MAC resulting from a specific choice of D from the MAC where all D from the 
feasible set are allowed, we denote the former as MAC(F[^1, Pi, F^^l, P2, G, D, Pr) and the latter as 
MAC(fW, Pi, F[2] , P2, G, {D} , Pr) 
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3.2.1 MAC [fW , , fPI , P2 , G, {D}, Pr) 

For the two user parallel AF relay MAC(fW,Pi,fP1,P2,G,{D},Pr), the received signals at the 
relays and the common destination are as follows: 

R = FWxi+FPlx2 + Ni^ (25) 
y = Tr[GD(|FWxi + FPlx2 + NR)]+n (26) 

Normalizing the noise to unit variance, the destination output can be expressed as 

Tr(GDFW) Tr(GDFPl) 

y = , ^ J= xx^ , ' X2 + n' (27) 

Vl + T>(D2G2) Vl + Tr(D2G2) ^ ' 

The power constraints are: 

E[x2] = Pi (28) 

E[xi] = P2 (29) 

Tr (^D^ (^/ + PiFW2 + P2F[212^) = Pr. (30) 

As in the point to point case, the power constraint can be substituted into the channel output 
equation to obtain a corresponding normalized form 

, _ Tr(GDFW)V^ 



y^TV [D2 (7 + PiF[i12 + P2FPI2 + PrG2)] 
TV (GDFPI) 

Tr [D2 (/ + PiF[i]2 + P2F[2]2 + P;^G2)] 



Xl 



X2 + n' (31) 



3.2.2 BC(G,Pr,fW,F[21,{D},P) 

We use the shorthand notation BC(G, Pr, F'^I, F[21 , D, P) to indicate the dual BC, i.e. broadcast 
channel with transmit power Pr, channel vector G from transmitter to relays, channel vectors F^^l 
and F[21 from the relays to receiver 1 and 2 respectively, relay amplification factor given by D and 
total transmit power used by the relays P. As for the MAC, we use {D} to indicate all feasible 
relay amplification factors are allowed and D to indicate a specific choice. For the dual broadcast 
channel, received signals at the relays and the two destinations are as follows: 

R = Gx + Nij (32) 
yi = Tr[FWD(Gx + NR)]+ni (33) 
y2 = TV[F[2lD(Gx + Ni?)] + n2 (34) 
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with power constraints 

E[x2] = Pn (35) 
Tr (I + PG^)) = P. (36) 

Notice that the relays are associated with power P and the source with power Pr. Normahzing 
the AWGN to unit power, 

Tr(FWDG) , , , 

y'l = I ^ ' =x + n', (37) 

'l + Tr (D2FW2) 

Tr(FPlDG) , ^ ^ 

2/2 = ^ ^ = re + (38) 

1 + Tr (DSFPla) 



Incorporating the relay power constraint into the normalized channel, 

Tr (fWdG) ^/P 

Vi = I x + n' 



Tr [D2 (/ + PFW2 + PrG2)] 

Tr(F[2]DG)/P , , , 

2/2 = ^ ' x + n' (39) 



TV [D2 (/ + PF[212 + P^G2)] 

Recall that for the point to point case, the dual channel and the original channel were identical. Even 
for the conventional MAC-BC duality, the channels on the MAC and BC are identical. However, the 
forms of the equivalent normalized channels for the MAC and BC presented above are significantly 
different. In particular, while in the BC, user I's channel is independent of F[2] and user 2's channel 
independent of in the MAC all channels depend on all parameters. A duality relationship 
between these two channels is therefore not a trivial observation based on their normalized channels. 
What makes the existence of a duality relationship even more unlikely is the fact that the channels 
themselves depend on how the total power is split among the users in the multiple access channel. 
With these apparent complexities, it is rather surprising that a precise duality relationship does 
exist between the MAC and BC described above, as we show in Section f5. II 

4 Capacity and Relay Optimization for Parallel AF Relay MAC 

Given a relay amplification vector D the capacity region of the resulting scalar Gaussian MAC is 
the well known pentagon that can be expressed as: 

C*^^^(FW,A,FPl,P2,G,D,P^)={(i?i,i?2) : <^^^{^ + P^T^^S)) ^ 
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R2 < log I + P2- ^ 



1 + TY(D2G2) 



R1 + R2 < 



' Pi [Tr (DGFW)]' + P2 [Tt (DGF[2])]' 
1 + TY(D2G2) 



Taking the union over all D that satisfy the relay sum power constraint gives us a characterization 
of the capacity region of MAC(F[i1, Pi, F^l, P2, G, {D}, Pr) as: 

CM^c^F[i]^P^^F[2]^P^^G^{D},Pk) = UDei.M.cC"'^''(FW,^i,F[^P2,G,D,<A) (40) 

where 

Vmac = {d : Tr (d2 (/ + PiFW2 ^ P2F[212)) = p^j (41) 

is the set of all relay amplification factors that satisfy the total transmit power constraint at the 
relays. Note that in the absence of any further characterization of D we are left with optimization 
over the entire space of feasible D, i.e. a R dimensional space. A brute force solution to such an 
optimization may be difficult as the number of relays increases. Theorem |21 solves this problem. 

4.1 Relay Optimization 

The following Theorem reveals the structure of the optimal relay amplification factor D for rate 
pairs on the boundary of the capacity region. 

Theorem 3 The optimal relay amplification matrix D to maximize any weighted sum of users' 
rates fiiR2 + ^2^2 with ^1, /U2 > 0, has the following form: 

D(e) = 7G (PiFW sine + P2F[2] cos^) (/ + PiFW2 ^ p^Y^2]2 ^ Pj^qA'^ (42) 



and 7 may be expressed explicitly as: 



Pr 



\^ Tr[G^ (PiFW sine + P2F[2l cose) (/ + P1FW2 + P2F[212 + P^G2) (/ + P1FW2 + p^F 



[2]2 



The proof of Theorem 13] is presented in Appendix lA. 11 Note that the optimization space over D is 
now only one dimensional, as opposed to the original R dimensional space. To identify a point on 
the boundary of the capacity region one only needs the corresponding 9. Therefore, the angle 9 in 
Theorem |31 plays a very important role. As we will establish in the following Theorems, |e| going 
from to tt/2 describes the boundary of the capacity region, with e = corresponding to the point 
where user 2 achieves his maximum rate, 9^^ (defined in Theorem corresponds to the points 
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where the sum rate is maximized, and |^| = 7r/2 corresponds to the point where user 1 achieves his 
maximum rate. This is also depicted in Fig|31 

Next we characterize the capacity region exphcitly by obtaining in closed form the rate pairs 
corresponding to the sum capacity as well as the individual user capacities. We will also present a 
system of equations whose solution is the rate pair (i?i,i?2) that maximizes fiiRi + fJ.2R2 for any 
positive /ii + /i2- We start with the extreme point /xi = l,/i2 = that characterizes the maximum 
possible rate of user 1. 

4.2 Maximum Individual Rate Ri 

Theorem 4 For the M^C(fW , Pi, F^l , P2, G, {D}, Pr), the maximum rate Ri that can be sup- 
ported is 

= log (^1 + PiPi?Tr ^G^fW^ J^/ + PiFW2 ^P2F[212 + PfjG2) "^^^ ^43^ 

This rate is achieved with relay amplification matrix D(0 = '/r/2), 

D(^/2) = 7P1GFW (1 + p^fW2 + P2F[2]2 + p^G^) (45) 



Further, if user 1 achieves his maximum rate Cio, i.e. D = D(7r/2), then the maximum rate that 
can he achieved by user 2 is 

1 + P2PR 



l^i=l 1 , o ^[112_Lp .[212^p„„2 ^^l^R Z^i = l 1^0 .[112^p .[2]2^ 



V ^^=1 i+p^^[112^p^^[212^^^^2 ' ^ - « ^^. = 1 i+p^^I112+p^^[212^p^^2 J J 

The proof of Theorem |1] is in Appendix IA.2I 

By symmetry, the expression for maximum rate 6*2^ for user 2 is obtained by switching indices 1 
and 2 in Theorem|lI Note that user 2 obtains his maximum rate when 9 = 0, i.e., with amplification 
matrix D(0). 

Theorem 13] leads to an interesting observation. Notice that C^^, i.e., the maximum rate for user 

1 in MAC(F[^1, Pi, F[2], P2, G, {D}, Pr) is not the same as the single user capacity of user 1 if user 

2 was not transmitting. In fact, for Pi,P2 > 0, 

^10 ^ c'^^^'(fW,Pi,G,{D},P^). (46) 

Recall that in the conventional Gaussian MAC (i.e. without AF relays) the maximum rate that 
user 1 can achieve is the same as his channel capacity as if user 2 is not transmitting. In the 
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conventional Gaussian MAC, even though user 2 transmits power P2, user 1 can achieve his single 
user capacity with power Pi if the rate R2 allows user 2 to be decoded and subtracted out before 
user 1 is decoded. The reason this is not true in the relay MAC is because of the power constraint at 
the relays. Even though user 2 can be decoded and his signal subtracted out at the destination, it 
does not make the system equivalent to one where user 2 was not transmitting at all. Transmission 
by user 2 affects the power constraint at the AF relays, thereby strictly reducing user I's channel 
capacity as compared to the case when user 2 is silent. 

Next we consider the sum rate capacity, i.e. /xi = /X2 = 1- 

4.3 Sum Rate Capacity 

Theorem 5 The sum rate capacity C^^ 0/ M^C(fW, Pi, fI^I, P2, G, {D}, Pr), i.e., a 2 user MAC 
with R parallel AF relays, when the user's transmit powers are Pi and P2, the total transmit power 
of all relays is Pr and the relay amplification factor is optimized over all feasible D : 

C^^ = log{l + SNR*) 

SNR* = ^ (^PiAn + P2A22 + [PiAu + P2A22)' - 4PiP2 (^11^122 - ^^2)) 

where 



^11 = V = Tr f G^fW^ (/ + p,fW2 + P2FP]2 + P^qA 

R 2f[2]2 / _,\ 

ttl + Pl/F + P2/P+PR5^ V ^ W 



R „2 .[1] .[2] 



^12 = V = Tr(G2FWFPl(/ + PiFW2 + p,FPl2 + P^G2)" 

^1 + Pi/WVP2/P + Ph52 V ^ ^ 

The sum rate capacity is achieved with relay amplification matrix Y){6 = 9^^), where 



"-tan-M ,.J^^±, ) (47) 



^SNR*-PrPiAu 
Further, the following rates pairs (Ri,R2) are sum rate optimal: 

«(2^1),pii(2^1)) = (lo, 



«(1^2),P^i(1^2)) = (^log(l + ^^^^^^l^),log(l + (l-/3)5#i2^) 

, . SNR* - P2PRA22 

where p 



2SNR* - PiPrAu - P2PRA22 ' 
and the rate pairs are achieved by successive decoding in the order 2 followed by 1, and 1 followed 
by 2, respectively. 
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Proof of Theorem [S] is presented in Appendix IA.3I Notice from Theorem [S] that the sum rate 
maximizing point is not unique. In other words the capacity region of the parallel AF relay 
MAC(fW,Pi,F[2],P2,G,{D},Pr) is not strictly convex in the region of sum rate optimality, at 
least for fixed Pi,P2. While the capacity region is not exactly a pentagon, it is similar to a pentagon 
with two smooth corners. The smooth corners arise from the different optimal D that maximize 
I11R1+ I12R2 for different ^1, /i2- This is similar in shape to the capacity region of a MIMO Gaussian 
MAC where different input covariance matrices are optimal for different /^i, /X2 pairs. Similar to the 
MIMO MAC, the capacity region is a straight line around the sum rate optimal region, connecting 
the two different rate pairs that both achieve the sum rate capacity. The intermediate points on the 
line correspond to rate pairs that are achieved with time sharing of the extreme points described 
in Theorem 121 

4.4 Optimal Rate Pair to maximize ^i^i + ^2^2 

Theorem 6 For any ^i,/i2 with > ^2, the weighted sum rate fiiRi + fJ.2R2 for the parallel AF 
relay M^C(fW , Pi, pl^l, P2, G, {D}, Pr) is maximized by the rate pair (Pf Pf given by 

Pf'^^ = log(l + 5iVP^), Pf''^^ = logfl + -^^') (48) 



1 + SNR*^ 

where SNR* and SNR2 are the solutions to the following simultaneous equations: 



IPiSNR"^ 



Proof of Theorem ini is provided in Appendix IA.4I Note that for /U2 > /Ui the corresponding rate 
pairs can be obtained by switching indices 1,2 in Theorem |SJ While the simultaneous equations 
do not appear to allow a closed form solution for general /Ui,/i2, closed form solutions can be 
obtained for the special cases of /Ui = l,/i2 = and ^ui = ^2 = 1 as discussed above, . In general. 
Theorem El allows a relatively straightforward numerical solution that is much easier than a brute 
force optimization over all feasible D. Solving the weighted sum of rate pairs optimization problem 
is especially useful for rate allocation problems as well as to plot the entire capacity region of the 
parallel AF relay MAC. 

4.5 The Capacity Region 

Combining the results of Theorems I3I4I5I and IBl the boundary of the capacity region for the parallel 
AF relay MAC(F[y, Pi, F^l, Pj, G, {D}, Pr) is graphically plotted by the following algorithm. The 
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R2 



A' 



A 9 = 




Ri 



e = 7r/2 
.F ^F' 



Figure 3: Capacity region of two user AF relay MAC(f[i], Pi, pl^l, P2, G, {D}, Pr) 

are labeled with reference to Fig. |51 

B: Draw a straight horizontal line from (Pi,i?2) = (0,C^^) to (i?i,i?2) = (C?\Cf ). 
C: For 6* : ^ 6^^ 



I cos ej (^I + PiFW^ + PaFl^]^ + P^G^) ^ 
[Tr(GDFW)]2 



D(6l) = G (^PiFWsin6l + P2FP] 

/ [Tr(GDFW)]^ , / r.A,, 

= + ^^^^ TV [D^ (I + P,fU + P.FI^I^ + P,G^)] + ^-^^[^ (^"^"^^ 

/ [Tr (GDF[21)]2 \ 

" log + P^PRr^ (J + p^F[il2 + P2FI212 + P^G2)] J 

l),pii(2-l). 



(52) 



D: Draw 



a straight line from (Pl^(l 2),P^i(l ^ 2) to (Pl^(2 
D^ E: For 9:6^^^ sgn(6l^i)f 

D(^) = G (^PiFW sin^ + PaF'^] cos^) + PiF^^ ^ p^F^^l^ + P^G^) "^ 



(53) 



Pi = log^l + PiPfl; 



[Tr(GDFW)]2 . 
^Tr [D2 (/ + P1FW2 + P2FI212 + P^G2)] ) 

( [Tr(GDFl2])]2 ^ . .^..^ 

^ + ^-^- Tr [D2 (/ + P1FM2 ^ ^ p,G2)] + ^-^-[^ (^"^'^0]^^ 
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• F: Draw a straight vertical line from (i?i,fi2) = (C|°,Cl°) to (-Ri,i?2) = (C|°,0). 

Fig. shows the typical shape of the capacity region. Optimal relay amplification factor is indicated 
on the figure in terms of the parameter \9\. 9 is equal to zero between points A and B, it changes 
from to 9^^ (as defined in Theorem 1^1) in the curved portion from points B to C. 9 is constant 
at 9^^ between points C to D. 9 changes from 9^^ to sgn(^^^)7r/2 as we traverse the boundary 
from E, and it is again constant at 9 = 7r/2 from the point E to point F. The coordinates of 
the points A, B, C, D, E, F are all known in closed form as given by the preceding results in this 
section. Point A' outside the capacity region indicates the maximum rate of user 2 with power P2 
if user 1 is not transmitting, i.e. Pi = 0. This is strictly higher than point A which corresponds 
to the maximum rate of user 2 with power P2 when user 1 is sending constant symbols {Ri = 0) 
that can be subtracted out at the receiver. As explained earlier, unlike the conventional Gaussian 
MAC, user 1 sending constant symbols with power Pi is not identical to user 1 silent, even though 
the constant symbols can be subtracted by the destination. This is because transmission by user 
1 affects the relay power constraint. Notice that AF relays are only allowed to scale the input, i.e. 
they cannot subtract the constant signal from user 1. Similarly, point F' indicates the maximum 
rate of user 1 if user 2 is not transmitting, which is higher than the maximum rate of user 1 when 
user 2 is transmitting (point F). 

5 Duality Relationships in Parallel AF Relay Networks 

In this section we examine the duality of the parallel AF relay MAC and BC. In the previous 
section we obtained the capacity region of the parallel AF MAC. Obtaining an equally explicit 
characterization for the parallel AF BC is much harder because, as we will see, without time- 
sharing among different D the BC capacity region is often non-convex. The notion of duality is 
particularly useful in such a scenario. Knowing the capacity region of the AF relay MAC would 
enable us to find the capacity region of the AF relay BC if we could obtain a duality relationship 
between the two. Establishing such a duality is the goal of this section. 

5.1 Duality of Parallel AF Relay MAC and BC 

The following two theorems establish the duality relationship between the parallel AF relay MAC 
and BC. 

Theorem 7 




(54) 
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Given any relay amplification matrix D that satisfies the power constraint on the parallel AF relay 
multiple access channel MA C{F^^\ Pi, F^'^\ P2,G,T), Pr), there exists a dual parallel A F relay broad- 
cast channel BC{G, Pr, fW , F^l , kD, Pi + P2) such that any rate pair (Pi, P2) that can be achieved 
on MylC(FW,Pi,F[2],P2,G,D,PK) can also be achieved on PC(G, Pr, FW, F^l , k{D}, P) . k is 
chosen to satisfy the relay sum power constraint on BC{G, Pr, F'^^ , F'^^ , kD, Pi + P2) . 

Theorem 8 

C^^(g,P«,fW,fP],D,p) =Up„p,>o,P,+P,=P C^^^^(fW,Pi,f[21,P2,G,^D,Pr) (55) 

Given any relay amplification matrix D that satisfies the power constraint on the parallel AF relay 
broadcast channel BG{G,Pr,F^^\F^'^\'D,P) and given a rate pair [Ri,R2) that is achievable on 
this parallel AF relay broadcast channel, there exist Pi,P2 > such that Pi + P2 = P and a dual 
multiple access channel M^C(F[-^], Pi, F^^l, P2, G, kD, Pr) such that the rate pair Ri, R2 is achiev- 
able on MAG{F^^\ Pi, F^^l, P2, G, kD, Pr). k, is chosen to satisfy the relay sum power constraint on 
MA C(F[il , Pi , fM , P2, G, D, Pr). 

The proof of Theorems Q and |H1 is presented in Appendix IA.5I 

Note that the duahty relationship is a strong duahty in the sense that the parallel AF relay 
MAC and BC are duals not only for the optimal relay amplification matrix D but also for any 
feasible D. This is an equally strong result as for the point to point case. Moreover, as in the 
point to point case, note that the powers of the relays and the transmitter are switched in the dual 
channel 

5.2 Capacity Region of the Parallel AF Relay BC 

The duality relationship described in the previous section provides a method to compute the ca- 
pacity region of the parallel AF broadcast channel. As shown in Fig. |l]the BC capacity region may 
be found simply as the union of the MAC capacity regions. Interestingly, while the MAC capacity 
region for fixed Pi,P2 was found to be convex in all numerical results, the BC capacity region is 
often a non-convex region. Note that the non-convex nature of the capacity region is because we 
have not considered time sharing between relay amplification factors. Indeed with time sharing 
between different relay factors, the capacity region can always be made convex. However, if time 
sharing between relay amplification factors is not allowed, the union of the MAC capacity regions 
for all different combinations of Pi +P2 = P, and hence the BC capacity region is often non-convex. 
While we have not shown the convexity of the MAC capacity region for any fixed Pi, P2, we expect 
it to be always convex. This hypothesis is supported by the fact that for fixed Pi,P2 the MAC 
sum capacity is maximized by a unique relay amplification D matrix, which indicates time sharing 
between different D may not be necessary to achieve any given point on the MAC capacity region. 
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Figure 4: BC capacity region as the union of MAC capacity regions 

5.3 Generalized AF MAC-BC Duality with Multiple Antenna Relays and Mul- 
tiple Hops 

The AF relay MAC-BC duahty shown in section I^TTl is not hmited to single antenna relays or two 
hop networks with parallel relays. The same duality relationship holds when some of the relays 
have multiple antennas. The duality also holds when more than 2 hops are considered. Thus, the 
MAC-BC duality holds for AF relay networks whether they are purely parallel (two hop), purely 
serial (multiple hops with a single relay at each hop) , or several parallel AF relay clusters connected 
in series to form a multihop relay network. It holds whether the relays are distributed with a single 
antenna at each relay or they are able to cooperate as as a multiple antenna node. The dual network 
in general is defined as the network where the MAC source nodes become the BC destination nodes 
and the MAC destination node becomes the BC source node. The power constraints are changed 
so that the link on each hop is subject to the same transmit power as in the dual network. Consider 
the two hop MAC. The source to relay channel (let us call it link 1), is subject to total transmit 
power Pi + P2 = P while the relay destination channel (link 2) is subject to total transmit power 
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Figure 5: 3- hop MAC-BC duality with distributed multiple antenna AF relays 

Pr. In the dual BC, the relay power becomes P so that the relay to destination channel (link 1) 
is still associated with power P. The source power on the dual BC is Pr so that link 2 is still 
associated with power Pr. In the general case, consider an m hop MAC where successive hops 
(starting with the source) are associated with sum power constraints Pi,P2, - - ,Pm respectively. 
Then in the dual BC, starting with the source, the power constraints will be Pm, Pm-i, • • • , -Pi so 
that Pi is the power of the last set of parallel relays. 

Since the proof technique is the same in all cases, we illustrate the proof with an example of a 
three hop relay network when multiple antennas may be present at some or all of the relay nodes. 
Multiple antenna AF relays can forward any linear transformation of their inputs. In other words, 
the amplification factor for a multiple antenna AF relay node can be any matrix with the same 
dimensions as the number of antennas at the relay. The overall amplification factor of each group of 
relays corresponding to a hop of the relay network can be represented then as block diagonal matrix. 
For the three hop example in this section, we use symbols A, B to indicate the relay amplification 
block matrices for the two relay hops, with the blocks representing the scaling matrices of multiple 
antenna relays. As a matter of notation, we denote the scaling matrices on the MAC as A^,Bj„ 
while we denote the scaling matrix on the BC as A5, B5. Let IrxI be a column vector containing 
all ones. Note that while we assume R is the number of relays at each hop for simplicity, the 
same proof applies even if the number of relays at each hop is different. We will use the following 
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definitions to simplify tlie notation for this section: 



G 




(56) 




= fWi^,xi 


(57) 






(58) 






(59) 



Thus, G,F^^^,Ff^',Nij 1 are column vectors. Symbol H indicates the MIMO channel between 
the two relay stages. Pri and Pr2 are the transmit powers of the two relay stages. The precise 
definitions of these symbols will be clear from Fig. Eland the following analysis. 
The output signals for the MAC and BC are: 

MAC y = G^Bm (uA^ (t'^xi + F^^'xa + N^j„i) + Nr^ + n (60) 



.2/1 = F™^A,(H^Bb(Gx + NH,i)+NH,2)+ni, 
y2 = F''^^A,(H^Bb(Gx + NH,i) +Nh,2) +n2. 



With noise power normalized to unity the MAC-BC are characterized as: 

MAC y = _ _ _ +n, E[xi\ = Pi,E[x2\ = P2 

'1 + g'^B^B^G + G^B^HA™AT;HrB^G 

Relay Power Constraints for MAC: 

Pi TV f A,„F™F™^A^) + P2TV f A^f'^f'^^'^A^) + Tr(A^A^) = Pr,. 



BC 



Pi TV (B™HA^Ff''Ff''^A^H^B^) + P2Tr (B^HA^Ff''F''''^A^H^B^ 

+ Tr (B^HA^ A^H^B^) + TV(B„B^) = Pr2 

F'^'^^AjjH-^BjG 
^/l+F'^'^AfcA^fF'^'+F'^'^AtHT^BtBlfHAT^F''' ^ 

Relay Power Constraints for BC: 
PffiTr (BfeGG^Bf ) + TV(B,Bn = P^ji 



PffiTr (a5H^B;,GG^B^HA^) + TV (Af,H^B;,B^HA^) + Tr(A;,Af ) = P = P^ + P2 
The following theorem states the duality result in this case. 

Theorem 9 For the 3-hop distributed multiple antenna AF relay MAC and BC described above 
and depicted in Fig. [31 

CM AC ^-p[i] ^ fPI, P2, H, A, Pri, G, B, Pr2) C C^^ (g, Pi?2, ^iB^, Pru^^ , F™, F^l, ^aA'^, Pi + 
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and 

C^^ (g, Pr2, KiB^, Pm, H^, fW , pt^l , ^2 A^, Pi + P2) 

= Up,,P,>o,Pi+P,=P C^^^(FW,Pi,F[2],P2,H,A,Pfli,G,B,Pfl2) (62) 

Notice that duality holds when = /tiA^ and = At2B^. The constants are determined by 
the relay power constraints. 

6 Conclusion 

We explored the capacity and duality aspects of AF relay networks. The MAC-BC duality known 
for conventional one hop Gaussian channel was found to be applicable to two hop communication 
over AF relay networks where some of the relays may have multiple antennas. A unique aspect of 
the AF relay MAC-BC duality is that the powers of the transmitter and the relays are switched 
in the dual network. With distributed single antenna relay nodes we determined the optimal 
relay scaling factors for the entire capacity region of the relay multiple access channel. Closed 
form expressions were found for the sum rate and individual maximum rates while simultaneous 
equations were found that can be solved to determine any rate pair on the boundary of the relay 
MAC. The capacity region of the relay BC was evaluated using duality as the union of the relay 
MAC capacity regions over different power splits between user 1 and user 2 while keeping the total 
power constant and equal to the total relay transmit power on the BC. 

We conclude with a word about the generality of the results. While we focus on real channel, 
signal, noise and amplification factors, the duality results and the associated proofs extend with 
almost no change to complex channels, signals and amplification factors as well. While the duality 
can be extended to multiple users as well, we have focused on the two user case to avoid the 
cumbersome notational aspects of considering multiple users in addition to multiple relays multiple 
hops and multiple antennas. The capacity characterization of the relay MAC for the complex case 
will also require phase optimizations of the relay scaling factors which may be non-trivial. Finally, 
we conjecture the duality results will hold for AF relay MIMO MAC-BC as well. However MIMO 
channels may require more work to establish equivalence of capacity regions as the structure of the 
transmitted signal covariance matrix on the dual MAC and BC will be different. 
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A Appendix 

A.l Proof of Theorem 131 



Consider the normalized MAC characterization 

y = ~^ xi H — ~^ X2 + n (63) 



/SNRi /SNR2 X2 ^ , 

where VSNRT = /Pr?%r7?^^^¥v and ^/SNRi = /^^%7?^^^ have been defined as such to 
ahow compact notation. 

Characterizing the capacity region of a two user MAC is equivalent to an optimization problem 
where a weighted sum of rates is maximized, i.e. max/xii?i + 112^2- Without loss of generality 
let us assume [i\ > ^2- For every choice of the relay amplification matrix D we obtain a scalar 
Gaussian MAC for which the capacity region is a pentagon and this rate pair is maximized for the 
corner point: 

(SNR \ 
1 + Y^sNk j ^^^^ 

Thus, we need to solve the following optimization problem: 

(SNR \ 
i + biN Ki / 

= max(/xi -/i2)log(l + SNRi) + ;U2log(l + SNRi + SNR2) (67) 
= max/x'ilog(l + SNRi)+;U2log(l + SNRi + SNR2) (68) 

where ^2 and ^'i = — IJ-2 are both positive. Thus, we need to maximize 

max(l + SNRi)''! (1 + SNRi + SNR2)''2 (gg) 

subject to the power constraint (|28() . 

We start with the Lagrangian formulation 

imx) = (i I P^iT.tid,g,fPr Y'W Pi(Etic^.g./f')^ + ^2(Eti4g./f)^ V^ 



A 



i:4(i+^i/r+P2/p)-p. 

.k=l 



(70) 



Setting the derivative ^'^g^'"^^ to zero yields the KKT characterization of optimal dj: 



(1 + SNR^)^'i(l + SNR^ + SNR^)'^2-i 



2^1 (Eti dkgkfHa^fr + ^MEti dkgkfnmn 
^ + Etidl9l 
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f[2]N2 



(1 + EtiW 

Absorbing the constants into A we have 

R R 



+ (1 + SNR^)^'i-^(l + SNR^ + SNR^)^2 
= \2di (1 + Pi/W' + P2/f' ') (71) 



/X2(1 + SNR^ 



PiiYl dkgkf^^)g^ff^ + P2{Y, dkgkft^)g^ff^ - (SNR^ + snr^)^^^, 



fe=i 



+ SNR^ + SNR^) 



R 



k=l 



k=l 



Xdk + + 



Ai2(l + SNR^ 



R 



R 



PiiYl dkgkfi^f + P2{Y, dk9kft^? - (SNR^ + SNR^) ^ gldl 



k=l 



+ /xi(l + SNR^ + SNR^) 



fc=i 

R 



R R 
PiiY.dkgkft'^f -^^KY.9kdl 



k=l 



k=l 



k=l 

= XPr 



/X2(l + SNR^)(SNR^ + SNR^) + + SNR^ + SNR^)SNRt = APr 
H2{1 + SNR^)(SNR^ + SNR^) + n[{l + SNR^ + SNR5)SNRt 



(72) 
(73) 



Pr 



(1 + P1/P + P2/P) 



+//2(1 + SNRt)(SNR^ + SNR^ 



//2(1 + SNR^ 



Fl Ft ~ 

PliYl dk9ki^)g^fl'^ + P2iY^ )<7./f 



k=l 



k=l 



R 



+ /xUl + SNR^ + SNR^)Pi(^ ')5i/f ' 



(74) 



fc=i 



/X2(l + SNR^)(SNRt + SNR|) + /.;(! + SNR^ SNR$)(SNRt) ^ ^^^[1,2 ^ ^^^[2]2 ^ ^^^2^ 



/X2(1 + SNR^) 



PliYl dkgkft%^fl'^ + P2iY^ )ft/f' 



k=l 



k=l 



+ 



R 



/.;(! + SNRt + SNR^)Pi(^ 45fe/i'^)5i/f ' 



(75) 



k=l 

U2(l + SNR^)(SNR^ + SNR^) + + SNR^ + SNR^)(SNR^) _ 

R M] 
K(l + SNR* + SNR*) + /.2(1 + SNR^)]Pi(j; ') ^^2 ^ [212 ; 

i+p,fr+P2fr+PRgf 



R 



+M2(l + SNRt)P2(J^45fc/; 



9if, 



2] 



fe=l 



m-2 m 



D = G (ciPiFW + C2P2FPI) (7 + PiFW2 ^ p^p[2]2 ^ p^^^) where ci , C2 > 0, 



(76) 



(77) 



24 



Defining 7 = y^c^~+~c| and tan0 = Ci/c2, the optimal D for maximizing + /U2i?2 can be 

expressed in the form: 

D = 7G (^PiPW sine + PaF'^l cos^) + PiF^^ ^ p^F^^l^ + PrG^)'^ . (79) 

Note that the constant 7 can be evaluated from the power constraint as in Theorem El ■ 

A. 2 Proof of Theorem m 

We are interested in the case ^1 = l,/i2 = 0. Substituting these values into the optimal D 
characterization of H76j) we have 

^« t^i l + A/™' + P2/f" + Pi?5f 

which corresponds to 9 = 11/2 and leads to 

R 2 f [i]2 
SNRt = PiPR^ (81) 

which gives C|*' = log(l + SNR^). In order to find C^^ , we need SNR2. From the definition, 

SNR^_P2(Eti45fc/P)' 



From (|8Uj) we have 



(82) 



(Eti45./f') fel + Pli^'^'VP2/fVP^,5.^ 

SNR^ = (I^ti^fcgfc/iy sNR^:^ (84) 



Ej=i 



SNR5 = ftFH^ ■^'•■/. 



Substituting into (|65() we have the expression for C. 
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A. 3 Proof of Theorem [H 

For sum rate we have ^11 = ^12 = 1, i-e. fi'i = 0,/X2 = 1. Substituting these values into (f7H|) and 
defining SNR'^ = SNR^ + SNR^, we have 

, SNR^ _ Pi (E grf/W) gfc/f ' + (E g^/'^O gfc/f .... 



and (5:5'i/[^0 ^ = PiA,, [Y^gdf^A + P2^22 (E^^^/'" 



Pb 
SNR* 

P^ 

Eliminating the factor gd/l^l) ^^"^^ ^^^^ equations leads to the quadratic equation: 

(Pl^n+P2^22)^^ + PlP2(AnA22-^?2) =0 (91) 



V Pr / -TR 



SNR* _ ^1^11 + ^2^22 ± V (Pl^ll + ^2^22)' - 4P1P2 (Aii^22 " Al^) 

^ "pT^ ^ ^^^^ 

The smaller root of the quadratic equation represents the minimum value of the objective function, 
while the larger root is the desired maximum value, which gives us the sum rate capacity expression 
of Theorem |SJ 

From definition, tan0 = |^ value is computed from H89() . 

Finally, we obtain the sum-rate optimal rate pairs by finding out SNR^ and SNR2 as follows. 

SNR* - A(Eg^/W)VP2(Egrf/t^Jf 

SNR - TTEdV ' ^ ^ 

Define (3 = — — ^ so SNR^ = /3SNR^ SNR^ = (1 - /3)SNR* (95) 
SNR 

Substituting the value of gd/W) f'^^™ (jS^ we have the definition of (3 as stated in Theorem [SjH 
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A. 4 Proof of Theorem IHl 

Continuing from ((7S)) . 

/i2(l + SNRt)(SNRt + SNR^) + + SNR^ + SNR5)(SNR^) 



di- 



Pr 

+ SNRi + SNR^) + ^2(1 + SNRt)]Pi(^ 4fffc/f ,12 m ; 

^ r 1 O f 

+/i2 (1 + SNRt ) P2 C£ '^kOkfT ) ma [212 o (9^) 

Multiplying both sides of the equation with and summing up all the equations for i = 1 to 

i = R, we have: 

/i2(l + SNR^)(SNRt + SNR^) + + SNR^ + SNR^)(SNR^) = 

^^2(1 + SNRt)(PRPiAn + PRP2Au/a) + fi[{l + SNRJ + SNR5)PRPiAn (98) 



^^"'^ " = ftlfi^ = V SiSfe- Similarly 



/i2(l + SNR^)(SNR^ + SNR^) + fi[{l + SNR^ + SNR^)(SNR^) = 

^2(1 + SNR\) {PRPiAua + PRP2A22) + + SNR{ + SNR5)P^jPiAi2a (99) 

Simple algebraic manipulation of these two simultaneous equations gives us the result of Theorem 

il ■ 

A. 5 Proof of Theorems [7| and HI 

Starting with the normalized BC as in H39() . without loss of generality we assume 
[Tr(FWDG)]2p [Tr (fP1DG)]^P 



Tr [D2 (/ + PFW2 + PrG2)] " Tt [D^ (/ + PFPIs + PrG^)] 

Thus, user 1 is the stronger user in this degraded AWGN broadcast channel. 

We wish to show that any rate pair {Ri,R2) achievable in MAC(fW, Pi, Ft^l, Pg, G, D, Pr) is 
also achievable in BC {G,Pr, FI^I , F^l , kD, Pi + P2) . For this we need the following lemma. 

Lemma 1 Consider the following rate pair {Ri,R2) on the boundary of the capacity region of 
MA C(F[i1 , Pi , fM , P2, G, D, Pr) with 

r^MAC ^ ,^J,, PiP^[rr(GDFW)]2 \ 

^ ^ Tr [D2 (/ + P1FW2 + p2F[2]2 + P^G2)] + P2Pfi[rr(GDF[2l)]2 J ^'^'^ 

P*^^^ - loJil P2PH[Tr(GDFl2])]2 \ 
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If this rate pair is achievable by any scalar AWGN broadcast channel where user 1 is stronger than 
user 2, then every achievable rate pair of MAC{Y^'^\ Pi,Y^'^\ P2,G,T), Pr) is also achievable in that 
broadcast channel. 



R2 



BC 



MAC 



Ri 



■45" \45° 



Figure 6: Proof of Theorem^ 

Proof: The proof of Lemma ^ is evident from Fig. |SJ The 45° slope hnes are tangent on a 
capacity region only at the sum rate maximizing point. For the MAC this is the slanted edge of 
the pentagon, while for the BC it is the point where all the rate and power is allocated to the 
stronger user, i.e. user 1 in the figure. As is evident from the figure, if the marked corner of the 
MAC capacity region is inside the capacity region of the BC then the entire MAC capacity region 
must be contained within the capacity region of the BC. ■ 
Therefore, we only need to show that the rate pair R2) specified in Theorem^is achievable 
in BC(G,Pr,FW,F[21,kD,Pi +^2)- Now, from the capacity of the scalar AWGN BC, we know 
every rate pair on the boundary of the capacity region of BC(G, Pr, fI^I, F'^I, kD, Pi + P2) can be 
written as: 

- l^' + Tr[D2 (/ + PFW2 + P^G2)]J ^^^^^ 

j,BC ^ . {l-a)PnP[Tr{Fi^^BG)]^ \ 

^ ^ 1^ aP^jP[Tr(F[2lDG)]2 + TV[D2(/ + PF[2]2 + P^G2)]^ ^ ^ 

Now, setting i?f = Rf^ we obtain the value of a 

PiTrrD2(/ + PFW2 + p^G2)l ^ ^ 

a = fl05) 

PTr [D2 (/ + P1FW2 + p2F[2]2 + P^G2)] + PP2Pr[Tt (F[21DG)]2 ^ ^ 

It is easily verified that < a < 1 so it is a feasible power split for the broadcast channel. 
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Similarly, setting iZf ^'^ = R^^ we again solve for a to obtain: 



(1 



PTr [D2 (/ + PiFW2 + p^Y\2]2 + p^Q^)] 



= PTr 



+ PTr 



a 



PiTr [d2 (/ + PFW2 + P^G2)] 



PTr [D2 (/ + PiF[i]2 + P2FI212 + P^G2)] + PP2Pi?[Tr (FPIdG)]^ 



(106) 

(107) 
(108) 



Thus, the same value of a satisfies both Rf^^ = Rf^ and R^^^^ = Pf^. Therefore, with 
this value of a, BC(G, Pr, FI^I, F^l, kD, Pi + P2) can achieve the rate pair (Pf Pf ^'^) spec- 
ified in Theorem ^ Applying the result of Theorem ^ this implies that the capacity region of 
MAC(fW,Pi,F[2],P2,G,D,Pr) is contained within the capacity region of its dual BC. This es- 
tablishes Theorem [7| 

Theorem |H1 is also easily established from the above result. This is because we have shown not 
only that the dual BC achieves the chosen rate pair on the boundary of the MAC, but also that 
this rate pair is on the boundary of the capacity region of the dual BC. In other words, we have 
shown above, that any choice of Pi,P2 on MAC(F[^1, Pi, f12] , P2, G, D, Pr) achieves a point on the 
boundary of the capacity region of BC (G, Pr, Fill , F[21 , kD, Pi + P2) . Moreover Pi , P2 are related 
to a by a continuous function. Now, as we take Pi from to P while keeping Pi + P2 = P, we 
trace the entire boundary of the capacity region of BC (G, Pr, F'^I , F[2] , kD, Pi + P2) . This implies 
the result of Theorem |H| ■ 



B Proof of Theorem El 

We begin by incorporating the power constraints into the channel definition. For the MAC, includ- 
ing the power constraint for the relays with total power Pfj2 into the channel we have: 



y 



yP^^B„HA,„F''^xi + /P^^B„HA„f''^X2 



B„HA™F™||2+P2||B„HA™F'^'|p+Pfl2||G^B,„HA™||2 + ||B.^HA„||2+Pfl2||G^B.^||2- 



+ n 



Next, incorporating the power constraint for the relays with total power Pri into the channel we 
have: 



y 



^^g^b^ha^f'^i + 



PriPriy^t ^[2] 

— G B^HAmF' ^X2 +n 

L\rn. 



(109) 
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where 

Am = PiPRillBmHAmF^ ^11^ + P2-PRi||BmHAmF' 'll^ + Pr2-Pri||G BmHAm| P + -Pri | IB^HA^ 

— T ||2|| — [I]ii2 — ii2ii — [2]ii2 — ii2ii ||2 

+-Pl^i?2||G Bmll ||AmF II +P2-PR2IIG Bmll ||Am,F II +Pr2||G Bmll ||Am|| 

+-Pl||Bm,||^||AmF' ^11^ + P2||Bm|P||AmF' '||^ + ||Bm|P||Am|P 

Similarly, incorporating the power constraint for the relays with power P into the broadcast 
channel we get: 

(/p)f'^J^a,h^b,g ^ 

yi = I _ _ _ x + ni, 

PK2||AbH^BbG||2 + ||AfeH^Bb||2 + HA^IP + P||fW^A,||2 + PUFW^AbH^B^IP 

and incorporating the power constraint for the relays with power Pri into the channel we get: 



yi = J^^F'^l'^AfeH^BfoGx + m, 



where 



aW = P«iP^2||A;,H^BfeG||2 + P^i||AbH^B,||2 + PP^i||F''l'^AbH'^B6||2 

+PPi^2||F™'^Ab||2||B,G||2 + P||F''l^Ab|p||Bfe||2 + P^2||Afe||2||BbG||2 + ||Ab||2||Bfe||2 
[2] 

and ?/2 and are similarly defined. 

Proceeding as in the proof of Theorem [7| let us assume, without loss of generality that user 1 
is the stronger user on the resulting BC, i.e., 

^||fW^a,h^b,g||^ > ^||f[^i^a,h^b,g|P (no) 

According to Lemma ^ we need to show that the following rate pair 

^MAC ^ logfl+ APmPi^2||G^B HAmF'V \ 
V A„ + P2Pm^'i?2||G B„HA„F^^||2y 

j^MAc ^ log (i + P2PmPR2\\G^B^UA^t'\\^ )^ ^^^^^ 



(113) 

is achievable and is on the boundary of the dual BC. 

Rate pairs on the boundary of the degraded dual BC with user 1 stronger than user 2 are given 

by: 

i?f^ = log(^l + «^^Jj^||FW^A,H^B,G|p^ (114) 

V Afl + aPP^iPffillF^'I'^AfeH^BbGlpy 
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« = =[2^ (116) 

It is easily verified that PA„ = PiA|,^' + P2Af'. Thus, the same value of a satisfies both Rf^^'-" = 
Rf-" and R2^^'~" = i?^*^ and the chosen boundary point of the MAC is also shown to be a boundary 
point of the dual BC. By the same arguments as the proof of Theorem [3 the duality relationship 
is established. ■ 
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